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We develop a DCA(PQMC) algorithm which employs the projective quantum Monte Carlo 
(PQMC) method for solving the equations of the dynamical cluster approximation (DCA) at zero 
temperature, and apply it for studying pair susceptibilities of the two-dimensional Hubbard-model 
with next-nearest neighbor hopping. In particular, we identify which pairing symmetry is dominant 
in the U-n parameter space {U: repulsive Coulomb interaction; n: electron density). We find that 
Px+y- {d^2_y2-) wave is dominant among triplet (singlet) pairings -at least for 0.3 < n < 0.8 and 
U < At. The crossover between d^2_j,2-wave and p^^+y-wave occurs around n ~ 0.4. 
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I. INTRODUCTION 

Since the discovery of high-temperature 
superconductivityA the purse of identifying the micro- 
scopic mechanism for unconventional superconductivity 
such as d-wave in cuprates or p-wave in ruthenates 
has been a main driving force in condensed matter 
physics. Even now, 19 years after the discovery of 
high-temperature superconductivity, it is an issue of 
hot debate whether the one-band (two-dimensional) 
Hubbard model, the simplest model for electronic 
correlations, is becoming d-wave or -further away from 
half-filling- p-wave superconducting at low temperatures. 

There are indications that there is indeed (i^2_j^2-wave 
superconductivity close to half-filling.^ For example, the 
functional renormalization group (fRG) predicts super- 
conductivity for weak couplingi^ Also numerical quantum 
Monte Carlo (QMC) simulations for finite size systems'* 
observe an enhancement of the d2;2_j,2-wave pairing. A 
promising method for addressing these questions is the 
dynamical cluster approximation (DCA) ^ an extension 
of dynamical mean field theory (DMFT))& which also 
takes non-local correlations and wave- vector dependences 
into account. This is essential for describing p- and 
d-wave superconductivity; DMFT itself only allows for 
s-wave superconductivity by construction. Such DCA 
calculations^ find dx2-y2-weLve superconductivity. How- 
ever, when solving the DCA equations numerically by 
conventional Hirsch-Fye QMC simulations^ one is re- 
stricted to rather high temperatures TAS Hence, a diffi- 
cult extrapolation to low T is necessary which is further 
hampered since (numerically) exact statements also re- 
quire an extrapolation cluster size Nc—^oo. 

Our paper will focus on a different parameter regime: 
Motivated by Sr2Ru04, many authors have recently ad- 
dressed the possibility of triplet superconductivity in 
the two-dimensional one-band Hubbard model with fi- 
nite next nearest neighbor hopping t' (the t-t' Hubbard 
model) at intermediate electron densities. The results 
have been controversial as regards the dominant pairing 
symmetry: fRG calculationsiSiiiii^ concluded that when 
t' > 0.3t—0At and the Fermi level is at the van Hove (vH) 
singularity, the system becomes ferromagnetic at suffi- 



ciently low T. While a d-wave phase spreads next to the 
ferromagnetic phase in the U-t' diagram for vH band fill- 
ings, a p-wave phase exists when going away from the vH 
band filling and for sufficiently large t' . Third-order per- 
turbation thcory^^ showed that triplet superconductivity 
is realized even when ferromagnetic spin fluctuation are 
not dominant. On the other hand, QMC,^^ FLEX,^?^ and 
DCAS calculations concluded that triplet pairing does 
not become dominant for intermediate filling even when 
t' is as large as OAt. 

It should be noted that fRG, third-order perturbation 
theory and the FLEX approximation are valid only for 
weak coupling; finite-size QMC is possible only for U <2t 
because of a serious negative sign problem^ and the 
DCA has been performed only for high Ts.^ Thus, the 
question concerning p- wave (triplet) superconductivity in 
the Hubbard model is not conclusive yet. 

In this paper, we introduce a new route to address 
this question, solving the DCA equations by an extended 
version of the projective QMC (PQMC) method.*^ This 
DCA(PQMC) approach mitigates the T-extrapolation 
problem of DCA(QMC). We concentrate on the crossover 
between d- and p-wave instability in the intermediate 
electron density range 0.3 < n < 0.8 and present results 
for the paramagnetic spectral function and the dominant 
pairing symmetry for U At and t' = OAt. On an equal 
footing with these pair susceptibilities, we also calculate 
the ferromagnetic and antiferromagnetic spin susceptibil- 
ity of the DCA cluster. 



II. t-t' HUBBARD MODEL 

The t-t' Hubbard model reads 

NN NNN 
H =-tY^ cl^Cja +t'J2 4aC]a +UJ2 '^^"a- (1) 

Here, and c^o- create and annihilate an electron with 
spin a on site i of the two-dimensional lattice; the first 
and second sum are restricted to nearest neighbors (NN) 
and next-nearest neighbors (NNN), respectively. In the 
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following, all energies are given in units of t, correspond- 
ing to a bandwidth of D = 8{t). 



for projection. Typically, we performed 10^ to 6 x 10^ 
QMC sweeps. To obtain G{iu}) from G{t), the maximum 
entropy method (MEM) is employed, as in Ref. 



III. DCA(PQMC) METHOD 

In the DCA,^ the Brillouin zone is divided into Nc 
patches, with a coarse-grained Green function G'(Kp,Li;) 
and self energy I]c(Kp, oj) for every patch p: 

G(Kp,c.) - [^o"'(J^P + - S,(K„c.)]-\ (2) 



Here, Gg is the non-interacting Green function and the 
k summation averages over all k-points of patch p which 
is centered around Kp; N is the total number of all k- 
points. G and Sc then define an effective cluster of An- 
derson impurities which can be described by the non- 
interacting Green function 



(3) 



or its Fourier transform ^Jo (^)xiX • This defines a cluster 
problem He given by 5o('^)xiXj and a local Coulomb 
interaction on every cluster site Xi. In the DCA, this 
cluster problem has to be solved self-consistently together 
with Eq. @. 

Here, we introduce a new cluster solver which is based 
on PQMC^^ and particularly constructed for zero tem- 
perature. Just as in Ref. 1^ T = expectation values of 
an arbitrary operator O are calculated as: 
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where Ht is an auxiliary Hamiltonian for which we take 
the cluster defined by ^o(^)xiX3 without Coulomb inter- 
action {U = Q). The local one-particle potential of Ht is 
adjusted self-consistently to yield the same n as the in- 
teracting cluster He- At least for A^c = 1, this yields the 
same asymptotic behavior of QoiT) and G(t) for large t. 
Hence, the influence of breaking time translational sym- 
metry at T = is reduced. This gives a smoother Green 
function in the vicinity of t = 9, but does not affect the 
results after long enough projection, i.e., sufficiently away 
from T = 9. 

As in Ref. ITsl the limit /3 — > 00 can be taken ana- 
lytically for this Ht- Then, the interacting Green func- 
tion G is obtained via the same updating equations for 
the auxiliary Hubbard-Stratonovich fields as for finite-T 
QMC. But the PQMC starting point is different: a T^O 
Green function with open boundary conditions defined 
for < T, r' < 9 = 16t,18t, discretized into L = 48, 64 
slices for A"c = 4x4 = 16 cluster sites. For the measure- 
ment of physical quantities, C — 8, 10, 12 central time 
slices are taken, and the remaining V time slices on the 
right and left side of the measuring interval are reserved 



IV. SPECTRAL PROPERTIES OF THE 
PARAMAGNETIC PHASE 

Now, let us turn to the DCA(PQMC) results. We 
performed calculations for n = 0.3, 0.4, 0.6, 0.8 and U — 
2t,3M,4:t. Throughout the study, we took t' = 0.4t and 
Nc — 4x4 =16 DCA patches. In the upper panel of Fig. ^ 
we present the one-particle spectral function A(K.p,uj) = 
— l/7rImG(Kp, w), i.e., the averaged spectrum over the 
six irreducible patches p of the Brillouin zone, as indi- 
cated in the inset. Patches 2, 3, and 4 have contribu- 
tions at the Fermi surface and hence a contribution at 
uj — 0. Thereby, ACK^^lu) has a particular sharp peak 
because the vH singularity is located at (tt, 0) and (0, tt). 
This sharp peak also reflects in the fully k-integrated 
A{uj) = J2p ^(Kp, w), shown in the bottom panel of Fig. 
n Let us note that the quantitative and qualitative fea- 
tures are very similar to those of DMFT (dashed line); 
there is no pseudo gap for t' = 0.4, n — 0.8, U = 4t in our 
r = DCA(PQMC) calculation. 



V. CALCULATION OF CLUSTER 
SUSCEPTIBILITIES 

To discuss the magnetic and pairing instabilities, we 
first calculate the corresponding susceptibilities 



X^P-(Q,ri-r2) = > e 



Xi ,X2 



(TrC^^|(ri)cx^^(ri)4^^(T2)cx^^(r2)), (5) 
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Xi ,X2 



5(Ki)5(K2) 



(TrCK,|(Tl)c_j^^^(Tl)c^K2i(^2)cJ^,|(T2)) (6) 

for the 4x4 cluster at self-consistency. Here, cko- = 
X^x '^Xo- exp(iKX) and g(K) is the form factor, i.e., 



.9(K) = 



1 for s-wave, 

\/2s\i\{Kx) for wave, 

\/2fi\i\[Kx + Ky) for Px+y wave, 

coii{Kx)^cos{Ky) for dx2-y2 wave. 



(7) 



Second, we calculate Imx(i-^) from x(t) by MEM and 
from this obtain the static susceptibility via the Kramers- 
Kronig relation 



x{t^ — 0)— / Imx(ti;)/ajdtj. 



(8) 



In DCA, the susceptibilities of the Hubbard model are 
calculated from the above cluster x by solving the Bethe- 
Salpeter equation;^ However, to this end it would be nec- 
essary to calculate susceptibilities for two r's instead of 
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FIG. 1: (Color online) Top panel: Spectral function A(Kp,Lj) 
for U — 4t,n — 0.8 and the 6 (out of 16) irreducible patches p in- 
dicated in the inset. The inset also shows the non-interacting 
Fermi surface as a line. Bottom panel: k-integrated spectral 
function A{uj) compared with that of DMFT. 



Eqs. © and 10 which would tremendously increase the 
numerical effort. Hence, we look at the cluster suscep- 
tibilities Eqs. |(SJ) and Q for simplicity. We can expect 
that this quantity already captures the essential features 
of the competition between ferromagnetic and antiferro- 
magnetic spin fluctuation, or that between superconduc- 
tivity with different symmetries. 



FIG. 2: (Color online) a) Susceptibilities for d^2_y2-wave 
pairing in r-space for U = 2t, n = 0.8 (crosses) and U = 4t,n = 
0.6 (circles). The arrows indicate different cutoffs Tc up to 
which xi^) information was taken into account in MEM for 
obtaining the lmx{(^) in part b) and c). Dotted, dashed, long 
dashed, and solid line are for the corresponding Tc-J, of the 
main figure. For U — 2t,n = 0.8 (panel b), the result is basi- 
cally Tc-independent; but for C/ = 4t, n = 0.6 (panel c) the peak 
at small ui is largely enhanced if Tc is increased. Consequently, 
the static susceptibility xi'^ = 0) is basically Tc-independent 
in panel d) and strongly Tc-dependent in panel e). In the fol- 
lowing, we indicate the problematic cases as in panel e) by an 
arrow. 



VI. INSTABILITIES OF THE PARAMAGNETIC 
PHASE 



Let us mention one more technical point: The decay 
of X in T-space can be underestimated because of a finite 
number of projection time slices V. In this case, the 
MEM calculation of Imx(w) depends on up to which Tc 
value the MEM input xi''') is considered, see Fig. 
Whenever we encounter a strong cutoff (rc-)dependence 
like in Fig. |2] (right inset) we take a short cutoff Tc (only 
four T points). With such a short cutoff Tc which is very 
close to the asymptotical behavior Tc — > 0, we certainly 
do not overestimate x(0)- However, if the weak decay 
of X in T-space is not due to the finite V but due to 
physics, we underestimate x(0) this way. In Figs. 13 and 
0] these more problematic data points, where we possibly 
underestimate x(0), are indicated by an arrow. 



Let us finally turn to the physical results, starting with 
the (inverse) ferromagnetic and antiferromagnetic sus- 
ceptibilities plotted in Fig. |3| as a function of T for three 
exemplary parameter values; finite-T data have been ob- 
tained by conventional DCA(PQMC). First of all, let us 
emphasize the obvious: DCA(PQMC) directly provides 
well for the T = susceptibilities without extrapolation 
-in very contrast to the finite-T DCA(QMC) for which 
it is difficult to forecast/extrapolate the low-T behavior. 
Turning to the physical results, we see that the ferromag- 
netic spin susceptibility becomes larger than the antifer- 
romagnetic one for n < 0.6. Note that the vH filling is 
at - 0.6 for t' = 0.4. In fact, the fRG studiesiii^ found 
a strong ferromagnetic instability for t' > 0.3 ^ 0.4 at 
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FIG. 3: (Color online) Inverse magnetic susceptibilities as 
a function of T for [/ = 2t,n = 0.3 (left), U = At,n ^ 0.6 
(middle), and U = 2t,n = 0.8 (right). The arrows indicate 
those data points where x is possibly underestimated (1/x 
overestimated) , as discussed in Fig. |21 
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FIG. 5: (Color online) U-n diagram showing the dominant 
pairing symmetry; for n = 0.4, U = 2t it is too close to call 
which instability is strongest. 
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FIG. 4: (Color online) Same as Fig.|3 but now for p^, Px+y, 
and d^2_y2-wave pairing as indicated in the right panel. 



vH band fillings. Thus, the qualitative tendency of our 
results and those of the fRG studies are consistent. 

Next, let us discuss superconductivity. Fig.0]shows the 
inverse static pair susceptibilities as a function of T. First 
of all, we find that d^^-y^-'wave (px+^-wave) is dominant 
among singlet (triplet) pairings for all {U, n)'s considered. 
For n = 0.8, d^2_y2-wave dominates over triplet pairings 
as expected. When we reduce n to 0.6, antiferromagnetic 
spin fluctuation become weak, see Fig. Nonetheless, 



the 



-wave instability is still larger than those for 



triplet pairings. This is consistent with finite tempera- 
ture DCA calculation of Ref. J for U = 3t,n = 0.67, and 
suggests c?2.2_j^2-wave ordering for low enough T -even 
for far away from half-filling. These findings are also in 
accord with QMC calculations for a finite-size Hubbard 
modelfi^ showing a dominant dx2-y2-wave instability for 
t' ~ 0.4, n - 0.67, and U = 2t. The fRG study Ref. [H 
reports also a ferromagnetic phase for t' = 0.4, n ~ 0.6, 
and U > t, while there is possibly a d-wave phase for 
U < t. On the other hand, third ordered perturbation 



calculations^^ have shown that triplet superconductivity 
becomes dominant even for weak coupling. We find such 
triplet pairing only for smaller n's, while the detailed 
structure of k-dependence of the vertex corrections may 
be essential which is possibly smeared out by the coarse- 
graining of DCA. 



VII. DISCUSSION AND OUTLOOK 

Fig. 13 summarizes our main finding, i.e., the dominant 
superconducting susceptibility on the cluster. We find 
that c?^2_y2-wave pairing is dominant in an unexpect- 
edly broad range of fillings n, whereas Px+y-^fave prevails 
only for small n. The transition between p^+y-'wave and 
dj.2_j^2-wave is at about rt~0.4. 

The t-t' Hubbard model with i'-0.4 and has 
been considered to be a good model for the 7-band of 
Sr2Ru04ii Our DCA(PQMC) results indicate however 
that there is no triplet superconductivity for these n and 
t' values, at least for the moderately strong values of U 
we studied and which are e.g. considered in local density 
approximation (LDA)-I-DMFT calculations>i2i 

Hence our results suggest that a nearest-neighbor 
Coulomb interaction V, as suggested in Refi^, or a multi- 
orbital model is necessary for an appropriate description 
of triplet pairing in Sr2Ru04 below 1 K. In contrast 
to the standard finite-T DCA(QMC), we can expect that 
DCA(PQMC) will open the door to such low tempera- 
tures -even for calculations with orbital realism. Simi- 
larly, we can hope that DCA(PQMC) will allow for more 
definite statements concerning superconductivity in the 
two-dimensional Hubbard model in the future. To this 
end, still a — > 00 extrapolation and the calculation of 
lattice susceptibilities is necessary. But at least the prob- 
lematic extrapolation to low temperatures is mitigated; 



PQMC is the direct route to T = 0. foundation (RA) and the Emmy Noether program of the 
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